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The phase-space volume of regions of regular or trapped motion, for bounded or scattering systems with two 
degrees of freedom respectively, displays universal properties. In particular, drastic reductions in the volume 
(gaps) are observed at specific values of a control parameter. Using the stability resonances we show that they, 
and not the mean-motion resonances, account for the position of these gaps. For more degrees of freedom, 
exciting these resonances divides the regions of trapped motion. For planetary rings, we demonstrate that this 
mechanism yields rings with multiple components. 
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The Cassini mission is providing unprecedented amounts of 
data on the structure of Saturn rings |1]. As usual, surprises 
and new questions emerge which call for theoretical under- 
standing. Beyond the intrinsic interest on planetary rings 
these are now the only examples of flat astrophysical systems 
for which such detailed data is available. This makes their 
analysis paradigmatic Jit], which serves as analog for other 
systems like galaxies and planet-forming disks. A first order 
understanding of the dynamics and processes in rings has been 
obtained, but important questions remain unanswered |20 . In 
particular, narrow planetary rings are eccentric and display 
multiple components (or strands), kinks, clumps and arcs (or 
patches); the F ring of Saturn is a magnificent case HI 12] • 
Therefore, these narrow rings pose interesting questions in 
dynamical astronomy and nonlinear dynamics related to their 
stability, confinement and structural properties. 

Recently it became clear that the phase-space volume asso- 
ciated to regular orbits for bounded systems with two degrees 
of freedom (DOF) displays universal properties as a function 
of a control parameter ]j3l |4|]. For scattering systems such 
quantity is related to the trapped orbits. In particular, there 
are locations where a drastic reduction of the phase-space vol- 
ume occupied by the regular or trapped orbits is observed, i.e., 
gaps. Here we show that these gaps are related to the occur- 
rence of stability resonances (SR), and not to mean-motion 
resonances (MMR). The MMR correspond to locations where 
the period of the particle is a rational of the orbital period 
of a satellite or planet. Examples include Kirkwood gaps in 
the asteroid main belt JH 0] and the Cassini division in Sat- 
urn rings 17], where very few particles are observed. Other 
cases of MMR exist where particles accumulate locally, as the 
Jupiter Trojans or the Hilda group ]8]. The SR are defined 
by a resonant condition in the (linear) stability exponents of a 
central stable periodic orbit. We show in this Letter that non- 
linear and higher-dimensional effects related to the SR yield 
rings with two or more components, or strands, as those ob- 
served in the F ring ][ll0]. We use the scattering approach to 
naiTow planetary rings ]9] and a disk rotating on a Kepler orbit 
as example ] 10]. Our results can be used in other cases, rang- 



ing from particle accelerators 111 111 to galactic dynamics ] 12 ] 



scattering dynamics of a (ring) particle in a planar system with 
some intrinsic rotation (external force). In the planetary case, 
this rotation comes from the orbital motion of one or sev- 
eral satellites. This rotation creates generically phase-space 
regions of dynamically trapped motion, in what otherwise is 
dominated by unbounded trajectories; "generic" implies that 
it holds for a wide class of rotating potentials lll411 . includ- 
ing gravitational interactions lls, 15]. For a uniformly circular 
rotation, the system has two degrees of freedom (DOF) and 
a constant of motion, the Jacobi integral. For these scatter- 
ing systems the organizing periodic orbits generically appear 
in pairs through saddle-center bifurcations, one is stable and 
the other unstable. The manifolds of the unstable orbit inter- 
sect, isolating a region where trapped motion is of nonzero 
measure if stable orbits exist. By changing the Jacobi inte- 
gral, the central linearly stable periodic orbit becomes unsta- 
ble and a period doubling cascade sets in. This scenario turns 
eventually the horseshoe (invariant set) into a hyperbolic one, 
thus destroying any region of trapped motion lll6ll . Then, in 
the Jacobi constant space, the regions of trapped motion are 
bounded. Particles with initial conditions outside these re- 
gions escape rapidly along scattering trajectories; those inside 
remain dynamically trapped. Consider now an ensemble of 
initial conditions of independent particles distributed over the 
extended phase space, which contains entirely one region of 
trapped motion. Then, due to the intrinsic rotation, the pattern 
formed by projecting the trapped particles into the X — Y plane 
at a given time forms a ring UM- The ring is typically narrow, 
sharp-edged and noncircular; this scenario is generic fl4]. For 
more than two DOF, e.g., a nonuniform rotation on a Kepler 
elliptic orbit, the ring displays further structure |9i]. 

The simplest example to illustrate the occurrence of rings 
is a planar scattering billiard: a disk on a Kepler orbit fllOll . 
This system consists of a point particle moving freely unless 
it collides with an impenetrable disk of radius d, which or- 
bits around the origin on a Kepler orbit of semi-major axis R 
(d < R). Collisions with the disk are treated as usual ifioll : no 
collisions lead to escape. For the disk on a circular orbit, the 
simple periodic orbits are the radial collision orbits. They and 
their linear stability are given by ^ 



where resonances and transport properties are important 
The scattering approach to narrow rings ]9] considers the 
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TtD3^j = 2+[{A(j)f{l-tan^ e)-4{l+A(j)tane)]R/d. (2) 

Equation ^ defines the value of the Jacobi integral J for the 
radial collision orbits in terms of 0, the outgoing angle of 
the particle's velocity after a collision with the disk. Here, 
= (2« — l)7r + 20 is the angular displacement of the disk 
between consecutive radial collisions and « = 0, 1 , 2, . . . is the 
number of full turns completed by the disk before the next 
collision. The period of the disk's orbit is = 2n (Oj = 1). 
Equation ^ provides the trace of the linearized matrix D.'S^j 
around the radial collision orbits. Being a two DOE system, 
periodic orbits are linearly stable iff |TrD^j| < 2. 

To understand the phase-space properties that define the dy- 
namics, in particular, for many DOE, we consider a relative 
measure of the phase-space volume occupied by the regions 
of trapped motion |l3l|4|]. This quantity is a function of some 
control parameter and tunes the horseshoe development l,16ll . 
Eor the disk on a circular orbit a good choice is J. However, 
this quantity is not conserved for nonzero eccentricity e, thus 
being useless for more than two DOE. A convenient quantity 
is the average time between consecutive collisions with the 
disk, (Af). The average is defined, for a given initial condi- 
tion (ring particle), over the successive consecutive collisions 
times; in addition, we consider an ensemble of them. In our 
numerical calculations we considered an orbit to be trapped if 
it displays more than 20000 collisions with the disk; the next 
200000 reflections were used to improve the statistics. 

In Eig. [Tfa) we show the histogram of (Af ) for the n = 
trapped region; Eig. [Tfb) shows a detail of the correspond- 
ing ring 191]. We note that Eig[lja) displays some structure, 
namely, some rather well localized positions where there is 
essentially no trapped motion in phase space or, at least, a 
drastic reduction of its size. This structure is similar to those 
in Ref. [3], and is universal Then, universality extends 

to scattering systems in terms of the phase-space volume of 
the regions of trapped motion. We note that the structure of 
Fig- Ea) also resembles the structure of the population his- 
togram of asteroids in terms of their semimajor axis [6], which 
uncovers the Kirkwood gaps and the role of MMR. 

Guided by the similarity with the asteroids, we ask whether 
the gaps in Eig.[Tla) are related to MMR. The question can be 
answered analytically for the radial collision periodic orbits 
which are the organizers of the dynamics. Their collision time 
is tco\ ~ A0 /(Bf/, hence the MMR condition is tt^ox/Tj = p/q, 
with p and q incommensurate integers. In Eig. [Tfa), we have 
indicated the location of the lowest mean-motion resonances 
in the (Af) interval of interest; the actual resonances are in- 
dicated in parentheses. We have also included the 29:44 res- 
onance because of its proximity to the position of the main 
gap. While these resonances are definitely not low-order res- 
onances, we note that some may be very close to the gaps, 
while others are not. Therefore, low-order mean-motion res- 
onances are not related to the reduction of the phase-space 
volume (gaps) of the regions of trapped motion. 

To understand the gaps of Eig.[TJa), we consider the phase- 
space structure defining a surface of section at the collisions 
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FIG. 1: (Color online) (a) Histogram of (At) showing the relative 
phase-space volume of the region of trapped motion for the scatter- 
ing billiard on a circular orbit. The continuous lines represent the 
position of the lowest MMR, indicated in parentheses, that occur in 
the interval of (Af); dash-dotted lines correspond to the lower-order 
SR. (b) Detail of the ring corresponding to this region of trapped 
motion. 



with the disk in a rotating frame (see [9]). Thus, the average 
return time to the surface of section is precisely (Af). Eig- 
ures|2]di splay the surfaces of section for three values of the Ja- 
cobi integral around the location of the main gap in Eig. \lla). 
Satellite period-three islands squeeze the central stable fixed 
point, reducing efficiently the region of trapped motion. At 
some value of J there is a bifurcation with the central point 
being unstable. After this value, the whole structure reap- 
pears with some symmetric inversion. This phenomenon is 
intrinsically nonlinear and is universal iITtIi . Note that the ex- 
istence of the period-three islands in the surface of section is 
not equivalent to a 1 :3 MMR. The period-three fixed points in- 
volve different collision times. Adding these individual times 
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FIG. 2: (Color online) Surfaces of section for the scattering billiard on a circular orbit for decreasing values of J. The sequence illustrates the 
reduction of the phase-space volume of the trapped region due to the 1:3 SR corresponding in Fig.lTJa), from left to right, to (At) ~ 4.1382, 
(At) ~ 4.1409, and (At) ~ 4.1445. The black/red curves are the stable/unstable manifolds of the unstable periodic orbit. The blue regions 
correspond to regions that represent more than 90% of the phase-space volume occupied by the trapped orbits. 



matches a multiple of the collision time of the central fixed 
point, which in general is not a rational of Tj. 

The gaps in Fig. [Tta) are thus related to a reduction of the 
phase-space volume occupied by trapped orbits due to satel- 
lite islands (see also [ISJ). In LI 7,1 it is shown that the passage 
through such 1 :3 resonances (Fig.|2]l leads universally to insta- 
bilities, while resonances of order higher than 4 do not induce 
instability; for the resonance of order 4 the behavior depends 
upon which contribution (resonant or nonresonant) dominates 
the normal form OlTlllSll . These results follow from the struc- 
ture of the resonant normal form, i.e., the relevant (nonlinear) 
contributions to the stability analysis of a central linearly sta- 
ble periodic orbit. We thus consider whether these resonances 
are related to the structure of Fig.lTta). We write the eigenval- 
ues of the linearized map D^j around the stable radial colli- 
sion orbit as A± = e*'". We are interested in values of a that 
are rational multiples of 2n, i.e., a/{2n) = p/q with p and q 
incommensurate integers. The eigenvalues X± are related to 
each other by complex conjugation; by consequence, the p : q 
resonance is related to the q — p : q resonance. Since these 
resonances involve the stability exponents, we refer to them 
as stability resonances. 

From the definition of the purely imaginary eigenvalues X±, 
\TtD^j\ < 2, we have cos a — TrD^j/2. This can be writ- 
ten in terms of the collision time fcoi = /(^d using Eq. (|2]l. 
In Fig. [TJa), we have indicated the lower-order SR as dash- 
dotted lines. The correspondence with the gaps is astonishing. 
Therefore, we attribute the gaps in Fig.dfa) to crossing a SR. 

The SR, and, in particular, the 1:3, have interesting con- 
sequences in the context of rings when we consider a small 
nonvanishing eccentricity of the disk's orbit. Note that, for 
nonzero £, the system is explicitly time dependent, has two- 
and-half DOF and cannot be handled with the usual tech- 
niques. In Fig. [3ja) we illustrate the phase-space volume of 
the trapped region for £ = 0.0001, indicating the location of 
some e = resonances as a guide, and in Fig. |3jb) we show 
a detail of the resulting ring. In comparison to Fig. \lla), the 
gaps at the 1:3 and 1:6 SR are wider, but the overall struc- 



ture remains. The histogram is divided in three distinct re- 
gions separated by those SR; £ tunes the width of the gaps. 
We construct the ring using different colors for the initial con- 
ditions in the different (At) regions. The ring displays two 
well-separated components. Each one is related to a different 
region in phase space; this is a consequence of the large gap 
opened by the 1:3 SR. These components are entangled and 
form a braided ring [9]. 

We note in Fig. |3jb) that there are only two independent 
ring components instead of three, as expected from the three 
disjoint intervals in (Af) defined by the 1:3 and 1:6 resonance. 
This follows from the fact that the 1 :6 resonance does not sep- 
arate enough the phase-space regions around it, so the pro- 
jections into the X —Y plane of these two intervals overlap 
and do not manifest two components. Yet, as observed in 
Fig. lUb), the coiTesponding ring component displays a clear 
separation of the ring particles, except for a thin common 
strip, depending on which side of the 1 :6 resonance they come 
from. For larger values of £ we have observed rings with 
three strands Ist]. Therefore, multiple ring components are 
obtained by exciting SR through nonzero eccentricity. Note 
that these properties follow from the higher dimensionality of 
phase space and nonlinear effects. 

The fact that the overall structure of Fig. [3a) is similar to 
the two DOF case and to the population histogram of asteroids 
may be an indication of universality for higher dimensions. 
Yet, our results are inconclusive for this issue: Fig.Oa) repre- 
sents the first results of the phase-space volume occupied by 
trapped orbits for a system of more than two DOF. 

To summarize, we have studied the phase-space volume oc- 
cupied by trapped orbits using a scattering billiard, a disk ro- 
tating on a Kepler orbit, and succeeded in relating some of its 
structure to the occurrence of SR. For two DOF the universal 
structure of the phase-space volume occupied by trapped or- 
bits is extended to scattering systems fslHl. For more DOF 
the qualitative similarity of Fig. ^a) and Fig.flja) is the first 
indication that universality may hold also for higher dimen- 
sions; yet, this issue remains open. Stability resonances are 
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FIG. 3: (Color online) Same as Fig. [T]for the billiard on an elliptic 
Kepler orbit with eccentricity e = 0.0001. The resonances indicated 
correspond to e = 0. In comparison to the circular case, the gap of 
the 1:3 SR is wider. This causes a true opening in the corresponding 
ring, thus forming a two-component ring. The black inner compo- 
nent of the ring corresponds to the region on the left of the 1 :6 SR; 
the gray (red) is associated to the region between the 1:6 and 1:3 res- 
onances. Note that the 1:6 resonance gap is not wide enough to have 
a projection with another (distinct) strand. 



defined by a resonant condition on the stability exponents of 
the linearized dynamics around a central stable periodic or- 
bit. These resonances, and not the mean-motion resonances, 
manifest locally as a reduction of the phase-space volume of 
the trapped trajectories due to nonlinear effects. While SR 
have a local manifestation, they have global consequences: 
Nonlinear and higher-dimensional effects (nonvanishing ec- 
centricity of the disk's motion) lead to an effective separation 
of the trapping region in phase space [Fig.[3]!a)]. In the con- 
text of planetary rings, if this separation is large enough, this 



yields a ring with two or more components or strands, which 
may entangle and form a braided ring. This provides a sim- 
ple explanation of recent observations of planetary rings with 
multiple components |2]. These results should be interesting 
beyond the context of planetary rings, in systems where reso- 
nances and the phase-space structure are significant; examples 
include particle accelerators 111 Itl and galactic dynamics il2ll . 
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